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Abstract 

We show that the growth of the size with the number of partons holds in a 
Thomas-Fermi analysis of the threshold bound state of DO-branes. Our results 
sharpen the evidence that for a fixed value of the eleven dimensional radius the 
partonic velocities can be made arbitrarily small as one approaches the large N 
limit. 
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1 Introduction 


In this paper we will study the semi-classical behavior of a nonrelativistic gas of DO- 
brane clusters. DO-branes are point-like solitons carrying the quantum numbers of the 
first massive Kaluza-Klein modes of the eleven dimensional supergravity multiplet. 
They are conjectured to be the fundamental degrees of freedom of M theory in the 
inhnite momentum frame and have a ten dimensional rest mass m = h/R, where 
R is the radius of the eleventh dimension. The large N limit of the supersymmetric 
quantum mechanics of DO-branes is as yet poorly understood. Our work is a hrst step 
in this direction. The DO-brane clusters are bound by the attractive potential arising 
from their two body interactions forming a threshold bound state we will refer to as 
the “D-atom”. The universal part of this attractive interaction is velocity dependent 
and its leading behavior is of the form 

For a large atom, the majority of the D-clusters have comparatively large quantum 
numbers so that we can use the semiclassical approximation to describe the D-atom. 
Since our WKB analysis only requires the universal interaction between D-clusters we 
will assume, with no loss of generality, that they carry no net spin but will allow for a 
possible degeneracy of rig DO-branes per cluster. Such a cluster is a D-parton carrying 
Hg/R units of eleven dimensional momentum. In what follows, we will estimate 
the growth of the size of the D-atom with the total number N of D-partons. The 
consistency of our analysis requires that the de Broglie wavelength of a parton with 
nonrelativistic mass m moving in the effective potential 0(r) of the remaining partons 
must vary only slightly over the characteristic scale of the system [^, i.e.. 



( 1 ) 


where A(r) = {27ih)/p{r) is the de Broglie wavelength, and r is the nine dimensional 
radial coordinate. The WKB approximation assumes ~ 0, and is thus invalid in 
the vicinity of the turning points of the effective potential. Let us scale distances by 
the eleven dimensional Planck length Ip, with r —> IpV. Expressed in these units, the 


WKB condition requires that for a potential of the form (p r ", we have r" ^ << 1. 


The spatial momentum of the D-atom is equipartitioned among its N constituents 
in the semiclassical approximation. Thus, the D-partons populate phase space uni¬ 
formly with one parton per phase space cell filling all available states up to some 
maximum momentum p{r). This procedure is independent of the statistics of the 
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D-partons. The total number of quantum states in a spatial volume element dV is 
given by 

• ( 2 ) 

where (327r^/945)|p|® is the volume of a spherical ball in 9d momentum space. 

The growth of the size of the threshold bound state with the number of D-partons 
has been estimated in a mean held approximation in |^. If the atom occupies a box 
of rescaled size L® with constant number density po, i.e., N = poL^, one hnds the 
scaling behavior: 

L ~ , (3) 

which is also in accordance with the holographic bound of one parton per spatial cell 
of Planck size area 0. We will show that the holographic growth of the size of the 
threshold bound state with the number of partons continues to hold when the partonic 
interactions are included, at least in the WKB approximation, thus sharpening the 
conjecture of Matrix theory. 

Let us use natural units setting h=c=l. The energy of any individual D-parton is 
E=^Rp^ — (j){r). Since this energy is assumed to be bounded from above we can write 
down a self-consistent relationship between the density and the potential within any 
spatial volume element of the D-atom: 

/’(’■) = . ( 4 ) 

where (pmax is the value of the potential at the boundary of the atom which we can 
set equal to zero. In a conventional 3d atom the electrostatic potential $ satishes the 
Poisson equation. Substituting for the electron number density, we get the Thomas- 
Fermi equation for the self-consistent potential in a large atom [|^: 

~ $3/2 . (5) 

In the case of the D-atom, the self-consistent potential is dominated by the veloc¬ 
ity dependent gravitational interactions of D-partons. We will use the analogous 
Thomas-Fermi equations for the D-atom to estimate the growth of its size with the 
number of D-partons. In contrast with conventional atoms where E ~ iV^/^ the nine 
dimensional D-atom displays holographic behavior: the D-parton energy density de¬ 
creases as the number of partons increases. 
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2 The effective potential 


We can infer the general form of the effective potential of D-partons from a dimen¬ 
sional analysis of the two-body interactions of DO-branes. Neglecting accelerations 
and time derivatives, the relative motion of a pair of nonrelativistic DO-branes in 
Minkowskian spacetime is obtained by expanding in powers of the held strength the 
DBI action i||5| [01171 : 


S = Tq J dt Tt det {rj^y + 27ia + fermions , 


( 6 ) 


where we have neglected commntator terms for a conhgnration of well-separated 
branes. The DO-brane tension is given by tq = [@, and we have nsed the nsnal 

relations between the Dp-brane tension Tp, the radins of the eleventh dimension R, 
and the dimensionless string conpling Qs'. 


R = gf%, Tp = {27r)-P{lp)-(P+^'>g, 


(p-2) 


(7) 


The one-loop approximation to this action is given by the snpersymmetric matrix 
qnantnm mechanics of SU{N) matrices, 0(f) = 


‘S'o = - ^10* X 


23' 


43 


( 8 ) 


redncing to the snpersymmetric dynamics of an ordinary coordinate world line, X{t) = 
{27ia')(j)rei{t), when the branes are well-separated as in the semiclassical regime. 

By a rescaling of the fermions, the action can be expressed as a loop expansion in 
the dimensionfnl parameter g = ^^ which plays the role of h in the semiclassical 

analysis. The leading terms in the expansion S = namely 

s = gldt + 0 (^ 7 °), 

determine the length dimensions of the helds: 

[X] = -1, [dt] = -1, [0] = -3/2, 

consistent with the snpersymmetry variations 

= *67*0, 6ifj = (X*7*)e . 


(9) 


( 10 ) 


( 11 ) 


It is amnsing that the velocity dependence of the nniversal interaction between 
DO-branes follows from this dimensional analysis when the leading term is a Coulomb 
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potential with a 1/r^ fall-off in nine spatial dimensions. Since [S'o] = —4 the velocity 
dependence must be of the form ~ We will set the coefficient in the partonic inter¬ 
action Hamiltonian equal to alp/R^, where a is an undetermined numerical coefficient 
proportional to the known value for ordinary DO-branes 0. Thus, 


0(r; v) 


alp |v(r)-v(r')d 

|r-r'p 


( 12 ) 


The semiclassical estimate for the growth of the size of the D-atom is independent of 
the precise value of this coefficient. 


3 The Thomas—Fermi scale 


We can now write down the nonrelativistic energy density describing the relative 
motion of a pair of D-partons with number density and velocity (p(r),v(r)) and 
(p(r'), v(r')). Integrating over the spatial volume gives the Hamiltonian: 

» = Ar)p(r)dV -‘41 p{r’)dVdV’ . (13) 

The dimensionful parameters in this Hamiltonian can be absorbed by the rescaling: 


Lr, V 




(14) 


This scaling renders energy dimensionless and can be motivated by the form of the 
space-time uncertainty relation in string theory with D-branes 0: 


AXAT>a 


(15) 


where a' = Ip/R- The velocity correlations are easily computed in the semiclassical 
regime. We have 


< v(r) >=< v(r') >= 0, <v2-(r)>=^n^ < (v(r) ■ v(r0)2 >= , 

(16) 

where n^(r) is the maximum allowed velocity in the volume element dV(r). In the 
semiclassical regime this velocity is self-consistently related to the number density 
through the phase space relation 

pir) = Cl^Mry . (17) 
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where C = ■ Performing the average over velocities we obtain 

< |v(r) — v(r')|"‘>= ^ + n^(r') + . (18) 

Substitution into the Hamiltonian yields velocity dependent nonlocal interactions. 
Such interactions can be rewritten in the form of a local held theory by introducing 
three auxiliary scalar helds, 0^, with equations of motion: 


A0i(r) -f u®(r) 

= 0 

(19) 

A02(r) ffw^^(r) 

= 0 

(20) 

A03(r) + u^®(r) 

= 0. 

(21) 


In simplifying these equations we have rescaled r — \J Varying the Hamilto¬ 
nian with respect to the velocity yields an algebraic equation relating the maximum 
kinetic energy density of the D-partons to the potential energy density stored in the 
three auxiliary potential helds, Ui = —(j)i{r)\ 

= 13 ^ 40 ^ + 11^202 ^ _ ( 22 ) 

As boundary conditions on the parton density we require that it be normalizable at 
spatial inhnity with a power law fall-oh faster than a Coulomb potential. Within 
the bulk of the atom, i.e., at the origin, the density can be assumed to approach a 
(positive) constant with small power law corrections. We have also found solutions 
to the diherential equations with a slow power law decay for which this constant van¬ 
ishes. The detailed analysis of the diherential equations together with the algebraic 
constraint equation is given in the appendix. 

It is easy to see that partonic interactions will induce corrections to the naive 
scaling behavior L ~ following from the relation p = po = NL~^. Consider 
the power law solutions ~ .^- 4 / 11 ^ ^ distances of order the size of 

the atom, r ~ L, we get an improved estimate for p, from which we infer the scaling 
L ~ Restoring the factors of Ip, the dimensionful size grows as 

L ~ , (23) 

for the growth of the size of the D-atom with the number of partons. Unlike con¬ 
ventional atoms where N plays the role of the nuclear charge Z, and E ZV®, the 
mean kinetic energy and number density of the D-atom decrease upon increasing the 
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number of D-partons. We can verify that the estimated atomic size lies within the 
WKB regime. For the above mentioned solution, the self-consistent potential 0 scales 
as At the scale of the atomic size, L, using eq.(23) we find the scaling behavior 

^ and the WKB bound is readily satished. 

The space-time uncertainty relation motivates another analogy drawn from atomic 
physics. Consider the Heisenberg relation for a highly excited single electron atom. 
The uncertainty AX scales as n^, the principal quantum number, so that AXAPr^nh. 
Let us derive an analogous result from the uncertainty relations in string theory. 
Assuming the scaling behavior L U X-2/81, with vAT = AX, from 

eq.(15): : 

AX AT >ll/R , (24) 

we get the scaling inequality: 

L^/v ~ {ly , (25) 

with N playing the role of the “principal quantum number”. 

4 Spin dependent forces 

We have neglected spin dependent forces in the Thomas-Fermi analysis. In this section 
we will show that in the semi-classical regime where the Thomas-Fermi approximation 
applies, all relevant spin-dependent terms in the two body interaction can indeed be 
ignored. 

As in sect. 2, if we rescale all helds so that there is an overall coupling g in front of 
the action, the one-loop term of the effective Lagrangian for the relative motion of two 
DO-branes will be independent of g and will have dimension —1. The supersymmetric 
partners of the interaction v‘"jr^ will be spin-dependent interactions at the one-loop 
level. In addition to the dimensional analysis, there is another rule that all the terms 
must observe. One may assign a quantum number 0 to X*, 1/2 to ip, —1/2 to e and 1 
to dt, consistent with supersymmetry and the form of the action. Thus the quantum 
number N = Nq + l/2Xj is the same for the same super-multiplet. This number is 
4 for n^/r’^. Together with dimensional analysis, this rule restricts the possible terms 
which can appear in the same supermultiplet as n^/r^. Schematically, we list the 
allowed terms as follows: 
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Now, the spin-flip does not occur without interactions, so that for the purposes of 
estimating the contribution of spin dependent interactions one can ignore terms con¬ 
taining derivatives of fermions. Thus, the relevant terms remaining from the above 
list are 

y3.02/^8, 

where the first term can be interpreted as a spin-orbital interaction, the second term 
as a dipole-dipole interaction, the third term as both, and the last term contains a pos¬ 
sible quadrapole-quadrapole interaction. The magnitude of both the spin-dependent 
terms and the spin-independent term can be expressed succinctly in the form 

(nr)"/r^^, where the spin factors can be dropped since they are independent of r at 
the characteristic scale of the large N bound state. Namely, {vL)^ ~ pL® ~ N, so 
that [vL]'^/ ~ It follows that this term grows with n, the largest value 

being n = 4, corresponding to the v^/r'^ interaction. Of course, we could improve 
upon this estimate by taking into account our previous Thomas-Fermi analysis. 

5 Conclusions 

We have shown that the holographic growth of the size of the threshold bound state 
with the number of partons 0 continues to hold for interacting D-partons, at least 
within the semiclassical approximation. We have thus sharpened the evidence for 
one of the key conjectures of Matrix theory, namely, that with R hxed as iV —> cx) 
the velocities can be made arbitrarily small at large N. We hnd that the qualitative 
behavior of the mean held analysis survives partonic interactions. This raises the 
interesting question of whether it is possible to compute the amplitude for the scat¬ 
tering of two D-atoms within the semiclassical regime. This would enable a direct 
test of longitudinal boost invariance within the semiclassical regime which, while it 
misses the tail of the distribution, is sensitive to the bulk of the kinematical phase 
space explored in the scattering. Boost invariance requires physics to only depend on 
the ratio of the sizes of the two bound states. It would be remarkable to understand 
how the holographic growth of particles can avoid contradicting longitudinal boost 
invariance in this picture. 
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6 Appendix 

We wish to hnd solutions to the nonlinear differential equation 

A0 = -// , (26) 


of the form (j){r) = Ar°'f{r) with S a positive rational number, where /(r) is regular 
either at the origin or at inhnity. Applying Frobenius’ method we can show that such 
solutions exist if and only if the leading power a takes the values 0, 2 — d, or — 
Substituting for the radial part of the d-dimensional Laplacian, = r^~'^dr{r'^~^dr), 
the function /(r) can be shown to satisfy the equation: 

+ (2a + d — l)r/'(r) J- a{a + d — 2)/(r) = —. (27) 

Let us expand /(r) in the Taylor series /(r) = 1 J- + • • •, where the in 

the limit r —0, and 6*^*^ < 0 in the opposing limit r — oo. With a{a + d — 2) = 0, and 
a=2a + d — l, we have a solution in the form of a recursion relation for the coefficients 
CA): 


=n{a{S-l) + 2), 


f,(»)(b(»)+a-l) 


2^m>l (5—m)!m! 


^n-m Q{ni) . . . (j{nm) 




( 28 ) 

where YAlLi Ui = n — 1. The boundary condition corresponds to the helds approaching 
a constant in the neighborhood of the origin, i.e., a = 0. At inhnity we can have a 
Goulomb potential with a = 2 —d. In addition, there exist solutions to the differential 
equation with a slow power law decay when a = —2/{5 — 2). In nine dimensions such 
a solution to the differential equation can have complex coefficients. 

The Thomas-Fermi equations for the D-atom take a similar form, with A^j = 
—/j 0 oS 00 = a positive dehnite function satisfying the constraint: 


- 1 ) + 9.).3 = 0 


( 29 ) 



Similarly, we can expand the velocity field in the neighborhood of the origin: 


00 — ^0 ( 1 — 


24 


4 ( 11 / 2 ) 

Ip 


9 11(1-A2)-26AiAo 


^2 _|_ 


( 30 ) 


with the relationship between coefficients of the auxiliary potentials l?>AiA^ + ll[A 2 — 
1)24o + 9^43 = 0. The velocity distribution can have a negative slope in the neighbor¬ 
hood of the origin. 
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